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Abstract 

In this paper we derive and present an application of hypergraphic 
oriented matroids for the purpose of enumerating the variable inter- 
dependencies that define the chemical complexes associated with the 
kinetics of non-linear dynamical system representations of chemical 
kinetic reaction flow networks. The derivation of a hypergraphic ori- 
ented matroid is obtained by defining a closure operator on families of 
n-subsets of signed multi-sets from which a "Z-module" is obtained. 
It has been observed that every instantiation of the closure operator 
on the signed multiset families define a matroid structure. It is then 
demonstrated that these structures generate a pair of dual matroids 
corresponding respectively to hyperspanning trees and hypercycles ob- 
tained from the corresponding directed hypergraphs. These structures 
are next systematically evaluated to obtain solution sets that satisfy 
systems of non-linear chemical kinetic reaction flow networks in the 
MAP Kinase cascade cell-signaling pathway. 



1 Introduction 

Numerous algebraic-geometric combinatorial models have been developed 
to analyze chemical reaction systems. Sellers [20] developed an algebraic- 
combinatorial model for identifying steady state and "cycle-free" solutions 
to discrete representations of derived dynamic systems, in which the reaction 
dynamics are given in terms of a linear transformation that defines set inter- 
section families of hyperplanes. The signed half-spaces implied by families 
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of sets of hyperplane arrangements, in turn, generate a signed convex poly- 
hedron that characterizes all possible flows modulo scalar multiples. Goss 
and Peccoud [9] demonstrated the use of Petri nets for studying chemical 
reaction systems. Oliveira et al. [12] showed that the Petri nets used by 
Goss and Peccoud are equivalent to directed hypergraphs (hyperdigraphs) 
and proposed using oriented matroids to analyze the topological structure of 
the hyperdigraph derived from a given chemical reaction system. Further, 
Oliveira et al. [13] E] demonstrated the utility of this topological hyper- 
graph analysis using the Krebs cycles, and then later in a sequel the EGFR 
signaling network. Yang et al. [22] also demonstrated the use of oriented 
matroids in analyzing steady state solutions to chemical reaction dynamics. 

The theory of oriented matroids provides a combinatorial algebraic- 
geometric abstraction of linear spaces based on signed orientations [U H] . In 
all generality, the stoichiometric number, which corresponds to the molecu- 
lar quantity of chemical reactant required to drive a reaction, will be greater 
than or equal to zero, and not restricted to the set {—1, 0, +1}. Therefore, 
we consider the stoichiometric number as being a "ranking" edge weight. 
This is the reason we have elected to define the oriented matroid represen- 
tation of the chemical reaction system over the Z-module which maintains 
the multiplicities. 

Multisets are used to track the multiplicity of member elements in a 
given base set S. Based on this notion, signed multisets are instantiated 
multisets with the added generalization that elements of the base set are 
allowed to have "negative multiplicity" . Multisets and signed multisets are 
used extensively in the study of hypergraphs and Petri nets [5] |T7] as a 
representation for vertex markings and edge weights. In the latter case, 
|17j uses the term "multirelation" , when the multiset base set is given as 
a relation S C D x D for some set D. We demonstrate that multisets 
are generalizable as a Z-module, over which we will construct a hyperdi- 
graph model of chemical reaction networks. An observable consequence of 
the corresponding matroid structure is a constructive method for obtaining 
"hyper-spanning" forests and enumerating the hypercycles in hyperdigraphs 
with weights [2| E] • This structure will enable us to qualitatively analysis 
of the steady state conditions in the chemical reaction network. 

Every hypergraph is a generalization of a graph [2 in which each hyper- 
edge may be incident with more than two vertices. The motivation for using 
hypergraphs as a representation of the chemical reaction network comes from 
the fact that not all of the reactions in such systems necessarily need to be 
unimolecular. A unimolecular reaction involves a single chemical species 
undergoing a molecular change to produce a different chemical species. In 



2 



general, many chemical species react as a chemical complex to produce a 
different chemical complex, which is in turn composed of many chemical 
species. We will therefore model the chemical reactions as directed hyper- 
edges, thus generating a hyperdigraph as a model of the chemical reaction 
network. 

We have validated and verified this discrete systems model approach 
by introducing the analysis of the mitogen-activated protein kinase cascade 
(MAP Kinase cascade) to demonstrate the effective use of applying the 
hyperdigraph model of chemical reaction networks [3]. The MAP Kinase 
cascade signaling pathway is a phosphorylation cascade that activates a set 
of cell regulatory molecules that govern cellular processes such as cell prolif- 
eration, differentiation, and development. It is interesting and worth noting 
from an evolutionary perspective that many distinct and disjoint cellular 
signaling processes have been identified that utilize the MAP Kinase cas- 
cade pattern. The MAP Kinase cascade is therefore considered to be an 
isoform for this category of cascading biochemical reaction networks |19j . 
The hyperdigraph representation provides a novel approach for the analysis 
of the steady state flux conditions. Further, we will utilize the hyperdigraph 
model of the MAP Kinase cascade to address an open question by jH] re- 
garding the graphic centrality rank index analysis of hyperdigraph models 
of chemical reaction systems. 

2 Hyperdigraph Model of Chemical Reaction Net- 
works 

First let S be a finite set of chemical species over which we will define 
chemical complexes. A chemical complex c G N s is a multiset defined on the 
set of chemical species, where the multiplicity c(s) (s G S) is the molecularity 
of species s in the given complex. We denote this finite collection of chemical 
complexes by C C N s . The finite set of reactions then is defined by the 
relation 1Z C C x C. Next define two projection functions p and tt from 
C x C into C, where p(c, d) \— ► c is to be referred to as the reactant and 
correspondingly, 7r(c, d) i— > d is the product of the reaction defined by the 
relation clZd . 

Now we define two matrices A, B £ N 7 ^ x,s ; where r = (c, d) S 1Z and 
s £ S, such that 

A(r,s) = p(c,c)(s) = c(s), 

and 

B(r,s) = ir(c,c)(s) = c(s). 
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Observe now that the matrices A and B correspond to the molecularities of 
the reactant and product chemical complexes, respectively, for each reaction. 
The net change in the species molecularity, of the given complex, can be 
modeled by the matrix N 6 'Z Sx ' R ' where 

N = {B- A) T . 

Further observe that for s € S and r G 1Z it follows that 

N(s, r) = B T (s, r) - A T (s, r) = 7r(r)(s) - p(r)(s). 

Therefore the columns of the matrix TV are presented as a finite collection 
of signed multisets defined on the set S. 

A hyperdigraph H = (V, £) is defined by a finite vertex set V and an 
edge set £ C {— 1,0, +1}^ such that each hyperedge E G £ is a signed set on 
the vertex set V so that E 1 partitions the set V into classes E + , E~ , and £7°. 
A weighted hyperdigraph H = (V,£, W) is a hyperdigraph TL = (V,£) with 
an assigned weight function given by W £ M^_ x£ . Evidently the incidence 
matrix of the hyperdigraph M £ M v ' x ^ is defined for v G V and E 6 £ by 
setting 

f ift;GE+, 
M(v,E) = E(v)W(v,E) = I -W(v,E), tiv€E~, 

0, otherwise. 

We observe that the multi-signed set model (S, 1Z) of a chemical reac- 
tion system defines a weighted hyperdigraph with incidence matrix N . The 
construction proceeds as follows. Let the chemical species of the chemi- 
cal reaction system be denoted by the finite vertex set S. We next obtain 
the hyperedge set £ by defining a mapping from the collection of reac- 
tions 1Z to { — 1, 0, +1} S given by £{r) = sgn(iV(-,r)) for r G 1Z so that 
£{TV) Q {-1,0,+1} S . Hence, each r in 7Z defines a directed hyperedge from 
which a definition for the weight function W 6 R + x ^ is given by setting 

W(s,r) = \N(s,r)\. 

From this realization is follows that the incidence matrix M G M s ' x7 ^ is 
defined by setting 

M(s, r) = £(r)(s)W(s, r) = sgn(iV(s, r))|JV(a, r)| = N(s, r). 

Hence, the incidence matrix N is shown to be sufficient to "model" the 
weighted hyperdigraph representation of the chemical reaction system de- 
fined by multi-signed sets on the set of chemical species. In this approach, 
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the directed hyperedges correspond to the set of reactions with the appro- 
priate associated weights, which denote the change in molecularities that 
are induced by the given reaction. 

Let /C G M + be the kinetic rate constants for the set of chemical reactions 
1Z. Also let X G be the concentration of the chemical species S. Now 
define a function p: 1Z — » M given by setting 

p(r) = n^) P(r)(s) , 

for every r in 1Z. The function p(-) is defined to be the chemical potential 
Next we define a function J: 1Z — > R by setting 



J(r) = K,(r)p(r), 

for all r in 7£ to define the flux. With these definitions we next define the 
linear operator • obtained by setting 

X = NJ, 

which corresponds to the differential equation derived from the dynamics 
of the chemical reaction network that was obtained from the constructed 
hyperdigraph model. 

The Michaelis-Menten ezymatic reaction will be considered for a 
small demonstration of the hyperdigraph model of chemical reaction net- 
works. The enzymatic reaction is kinetically represented as 

s + e , c p + e, 

where reaction r\ is complex formation of the substrate s with the enzyme e 
to produce the complex c; reaction ri is the disassociation of complex c into 
the substrate and enzyme, which is also considered to be the reverse reaction 
of r\ ; and reaction r% is the molecular reaction wherein the enzyme e disasso- 
ciates from the complex c producing the product p. The chemical species set 
is S = {s, e, c, p} over which the collection of chemical complexes is defined as 
C = {(1, 1, 0, 0), (0, 0, 1, 0), (0, 1,0, 1)}. We further define the relation K rep- 
resenting the three reactions by (1, 1, 0, 0)K(0, 0, 1, 0), (0, 0, 1, 0)^(1, 1,0, 0), 
and (0, 0, 1, 0)7£(0, 1, 0, 1). Finally, the incidence matrix is given by setting 

H r 2 r 3 
s /-ll \ 
e -111 
c 1-1-1 
p \ 1 J 



N 



5 



The hyperdigraph for this chemical reaction network is presented in Figure 



1(a) the vertices (chemical species) are denoted by circles and the directed 
hyperedges (chemical reactions) are denoted by rectangles. The arrows de- 
note the directionality of the hyperedge with respect each of its incident 
vertices. Thus, the hyperedge labeled rl represents transport of molecular 
mass from vertices s and e to vertex c. 




(a) Enzymatic chemical reaction. (b) A hypothetical example of 

a hyperdigraph. 



Figure 1: Two examples of hyper digraphs. The solid arrows indicate a 
hyperspanning tree. 



We now construct an oriented matroid, which we refer to as a hyper- 
graphic oriented matroid, on the hyperdigraph and show that the hyper- 
graphic oriented matroid is isomorphic to the steady state solutions of the 
chemical reaction network. It is worth noting that the hypergraphic oriented 
matroid corresponding to the steady state solutions of the dynamical system 
are equivalent to the ^-invariants of the Petri net chemical reaction model 
[9]. This equivalence is due to the fact that hyperdigraphs are isomorphic 
to Petri nets [T3] . 



3 Hypergraphic Oriented Matroids 

Let S be a finite set and consider the multisets x G Z s defined on S, where 
x(s) is the multiplicity of the element s in S. Denote to be the multiset 
that maps every element s to 0. We observe that allowing negative multi- 
plicities is a generalization of multisets [T7] , which are traditionally defined 
by mappings to the natural numbers [3 [6|. We will define the operation 
of addition on X s where, for arbitrary multisets x and y in Z 5 and s in 
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S we have (x + y)(s) = x(s) + y(s). Further, we define the unary opera- 
tion of negation as (— x)(s) = —x(s). It is now clear that defines 
an additive abelian group. Further we may define scalar multiplication as 
the function Z x Z 5 i— ► Z 5 , given by (ax)(s) = ax(s) for arbitrary a in Z, 
x in 7j S , and s in S. We observe that for a and /? in Z, and x and y in 
Z 5 , we have the following distributive properties (a + /3)x = ax + /3x and 
a(x + y) = ax + ay. We may therefore consider Z 5 as a module which is 
referred to as the Z-module [THl EE] . 

Recall the multiset model of the chemical reaction system defined pre- 
viously, wherein the set S is the set of chemical species and the matrix N 
represented a finite collection of multisets N(-,1Z) on S. Given that the set 
1Z corresponds to the hyperedges of the hyperdigraph, we have that the two 
dual oriented matroids correspond, respectively, to the linearly independent 
and linearly dependent subsets of 1Z. The minimal linearly dependent sub- 
sets of TZ will correspond to hypercycles of the hyperdigraph. Similarly, the 
maximal linearly independent subsets of 1Z will correspond to the hyper- 
spanning forests of the hyperdigraph. 

We will now proceed to construct the hypergraphic oriented matroid 
on the Z-module M C Z^. An operator cl is on the Z-module M as the 
function from 2 M into 2 M given by, for any X C M, 



It was shown by [7] that the function cl is a closure operator and so satisfies 
the following properties |15| . 



(CL1) iflCM, then X C cl(X); 

(CL2) iflCFCM, then cl(X) C cl(F); 

(CL3) iflCM, then cl(cl(X)) = cl(X); 

(CL4) if X C M, x G M, and y G cl(X Ux)\ cl(X), then x G cl(X U y). 
We observe that for X subset of M, cl(X) is a submodule of the Z-module 



For the construction of the matroid bases it is important to define an 
irreducible multiset. To do so, we will define a function g-.M^TL given by 
setting, for any x in M, 




m G M bm=Y^ a(x)x, for a G Z x , / b G Z 



M 0. 



g(x) = a(r)x(r) 
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where a in Z 7 ^ is chosen such that g(x) = gcd{x(r) | r E TZ}. We then define 
the mapping °: M — > M, which will be referred to as the "reducing map" by 
setting 

x(r) = ±- )X {t), 

for every r in TZ. A multiset x is irreducible if x = x. The intention of the 
reducing map is to obtain an analog of mapping a vector in a vector field to 
its corresponding unitary vector. Define to be irreducible. For a subset of 
X of M, we define X = {x E M \ x E X}. 

We now present an algorithm similar to Fourier-Motzkin [16j to construct 
bases for the dual hypergraphic oriented matriods M and M* of Z 7 ^. 

Algorithm 1. Basis: Let S be a finite set and let 1Z be a finite subset of 
li S , which will be represented by the integral matrix N E Z 5 * 7 ^. We will 
construct the mappings F £ Z 7ex ( s07e ) and F* £ Z 5x ^° 5 ) by 

s n n s 

F = k( N t Id(K) ) , F* = S ( N Id(S) ) , 

where Id(lZ) and Id(5) are identity matrices. We will then perform Gaus- 
sian elimination on F and F* where the row elimination process of using 
F(ri, ■) to eliminate F(rj,-) is 



c 
a 
b 



lcm{F(ri,s),F(rj,s)} 
c/F(ri, s) 
c/F(rj,s) 
bF{r J ,-)-aF{r ir ) 



After the completion of Gaussian elimination on the matrices F and F* , 
each is partitioned into submatrices as 



S 



S K 

f im(iV) Dom/ker(iV) \ 

v ker(iV) J 

TZ S 

im(iV T ) Dom/ker(7V T ) 

ker(iV T ) 



The respective basis sets for M and M* are then defined by the submatrix 
of F labeled ker(N) and the submatrix of F* labeled im(N T ). That is, 

B = {F(r, TZ) | F(r, S) = 0, re1Z}, 
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and 

B* = {F*(s,TZ) | F*(s,K) ^ 0, s G 5}, 

W^e i/ms obtain the sets M = (cl(-B)) and M* = (cl(B*)), which correspond 
to the desired hypergraph oriented matroids. 

We remark that even though B may not be a basis set for the oriented 
matroid, it is sufficient to identify a fundamental set of independent hypercy- 
cles as well as the matroid rank dimension. As such, future considerations 
may be to determine whether B is in fact a basis for M and also check 
whether X C cl(X). 

4 Hypercycles and Hyperspanning Trees 

The dual matroids M and M* correspond respectively to the cycle space 
and co-cycle space of the chemical reaction network represented by the in- 
cidence matrix N. The elements of M are hypercycles, the basis elements 
of which are minimal linearly dependent subsets of the reactions that define 
the hyperedges of the network. We use the term "hypercycle" to make a 
distinction between the elements of M and cycles in hypergraphs as defined 
by Berge [SJ. 

In the general theory of hypergraphs [2] , a chain of length q in a hyper- 
graph is a sequence v\ ri ► V2 rz * • . . . — q —* v q+ \ such that 

(CI) vi, V2, ■ ■ ■ , Vg are all distinct vertices of H; 

(C2) n, T2, ■ ■ ■ , r q are all distinct hyperedges of TL; 

(C3) v k ,v k+ i G supp(r fc ) with p(r k )(v k )p(r k )(v k+1 ) = = TT(r k )(v k )ir(r k )(v k+1 ) 
for k = 1, . . . , q. 

If q > 1 and v q +i = v±, then this chain is called a cycle of length q [2J. This 
general definition of cycle corresponds to the cycles which exist in the bi- 
partite representation of the hypergraph, and therefore does not incorporate 
the nonlinear dependencies that define the hyperedges in chemical reaction 
networks. For this reason we will use the term "closed loop" in leu of "cycle" 
for this class of hyper digraphs and define the more restrictive hypercycle. 

A hypercycle of length q is a multiset y G Z 7 ^ such that N(s, •) J_ y for 
every s in S with |supp(y)| = q. Observe that y(r) is the resulting mass flux 
caused by reaction r in TZ. Further, there is not an order restriction on the 
set of hyperedges that construct the hypercycle. Since N(s, •) • y = is the 
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net change in mass for chemical species s resulting from y, the hypercycles, 
and thus the matroid M, represent steady state flows in the network. 

A subhyperdigraph of a hyperdigraph TL = (V, TV) is a hyperdigraph 7i' = 
(V, 1Z') such that V and TZ' are subsets of V and 1Z, respectively, and each 
hyperedge in 1Z 1 is a multiset on V with N(V \ V , r) ^ for every r £ 1Z'. 
A hyperspanning forest in a hyperdigraph 7i is a maximal subhyperdigraph 

= (V,7?/) such that Ti.' does not contain any hypercycles and V = V . 

Corollary 1. The hypercyclomatic number for a hypergraph TC = (V,7Z) is 

c = rank(M) = nullity(iV) = \K\ - rank(iV T ), 
where M and M* are the matroids for ker(iV) and im(N T ), respectively. 
In general there may exist closed loops in hyperspanning trees as indi- 



cated in Figure 1 (b) , where the hyperspanning tree is depicted by the solid 
arrows. There are three independent closed loops in the hyperspanning tree, 
namely 







V5 




Vl, 




V3 


r-i 


V2 










V2 


ri 


V 3 


r3 


V2- 



V 5 V 3 V 2 " V 5 , 



The only hypercycle in this example is y := (0,0,1,1,1). The incidence 
matrix N is 



Tl T2 T3 T4 T5 



Vl 


( 


1 


-1 











\ 


V2 







-1 


1 


-1 







V3 










-1 


1 







V4 













1 


-1 




V5 


V 


-1 


1 





-1 


1 


/ 



5 MAP Kinase Cascade 

The MAP Kinase cascade is composed of six coupled enzymatic reactions 
and a complex formation [3], which are listed in the following tableau, in 
which in the interest of space we have abbreviated the enzymatic reactions 
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S + E 



C 



E + P asS 



P. 



PKC 

Raf : Raf* 
PP2-A 

GTP.Ras.Raf* 
MAPKK " MAPKK* 



MAPK 



PP2-A 
MAPKK** 



± MAPK-, 



tyr* 



MKP1 

Raf* + GTP.Ras " GTP.Ras.Raf* 



MAPK* 

Raf* : Raf** 

PP2-A 

GTP.Ras.Raf* 

MAPKK* s MAPKK* 



MAPK 



PP2-A 
MAPKK** 



tyr* 



MKP1 



± MAPK tyr * 



The hyperdigraph representation of the MAP Kinase cascade is presented 

in Figure [2j with the solid lines indicating one of many possible the hy- 

perspanning trees. The rank of the hypercycle, hypergraphic oriented ma- 

troid is 19, where there are three irreducible hypercycles for each of the 

Ei ^ 

six coupled enzymatic reactions i) S , P, ii) S + E\ S:E\, and 

E 2 

iii) P + E2 s P: E2; and one irreducible hypercycle for the reaction 

Raf* + GTP.Ras ; ^ GTP.Ras.Raf*. Recall that the hypercycles cor- 
respond to the steady state solutions of the dynamic system. As such, 
any linear combination of the 19 hypercycles is also a steady state solu- 
tion. The dual hypergraphic oriented matroid, which corresponds to the 
hyperspanning trees, also has rank 19. We observe that there are 8 closed 
loops within the hyperspanning tree presented in Figure [2] Six of the eight 
closed loops are linearly independent and are of length 2. The remaining 
two closed loops are each of length 4 and are linearly dependent on the six 
closed loops. Each of the six linearly independent closed loops correspond to 
one-half of a coupled enzymatic reaction, wherein the closed loop explicitly 
denotes the conservation of the enzyme while producing the product from 

the substrate as in the general depiction S + E «- S: E E + P. 

The two linearly dependent closed loops of length 4 correspond to a con- 
servation of the enzyme in a sequence of enzymatic reactions in which both 
reactions are catalyzed by the same enzyme as in the general depiction 

S + E S: E >■ E + S S : E P + E. 

Oliveira et al. |12| [13] [14] proposed and successfully implemented a 
method for rank ordering the chemical species and reactions within a chem- 
ical reaction system by counting the occurrences of each species (reaction) 
within the complete list of closed loops (Berge cycles) contained in the hyper- 
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Figure 2: Hyperdigraph of the MAP Kinase Cascade. The solid arrows 
indicate one of the many possible hyperspanning trees. 
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digraph of the chemical reaction network. It is observed that rank ordering 
the chemical species based on the complete listing of closed loops is anal- 
ogous to a rank ordering based on the centrality of the chemical species 
within the hyperdigraph [8j. Estrada and Rodriguez- Velaquez JB| character- 
izes the centrality rank index based on the hypergraph adjacency matrix. 
Using the hyperdigraph model proposed in this paper we can construct the 
hyperdigraph adjacency matrix for the chemical reaction network by setting 
L = A T B, where the matrices A and B were defined previously. 

There are 1456 irreducible closed loops within the MAP Kinase cascade 
hyperdigraph shown in Figure [2} Figure [3] shows the proportion of the 1456 
irreducible closed loops incident with each of the chemical species. For the 
chemical species the mean proportion is 0.538 and the standard deviation 
is 0.289. Oliveira et al. futher demonstrated that chemical species which 
occur most and least often are critical to the functionality of the chemical 
reaction network. Those species that occur most often within the closed 
loops, that is being most central, may be considered critical to the network 
as a "pinch point" or being utilized often in the signaling pathway. Those 
species that occur least often within the closed loops, that is being least 
central, may be considered critical to the network as an initiator or trigger 
for the cell-signaling process. 

For this analysis the designations of high and low centrality will respec- 
tively be considered by the values of the mean + standard deviation (0.827) 
and the mean - standard deviation (0.249). The most central chemical 
species are MAPKK** (0.926), Raf* : MAPK* (0.915), MAPK* (0.905), 
Raf* (0.885), and PP2 - A (0.856). The least central chemical species 
are GTP.Ras (6.87 x 10~ 4 ), PKC (1.37 x 10" 3 ), Raf (0.114), Raf:PKC 
(0.115), and Raf** (0.207). 

6 Discussion 

We have effectively demonstrated the construction of a multiset model for 
chemical reaction networks, which canonically generates a hyperdigraph. 
The hyperdigraph model was utilized in the qualitative steady state anal- 
ysis of the Map Kinase cascade, wherein it was shown that each of the six 
coupled enzymatic reactions contribute three independent hypercycles with 
an additional hypercycle corresponding to the complex formation Raf* + 

GTP.Ras N GTP.Ras. Raf* . The hypercycles are invariants in the dy- 
namical system of the chemical reaction network. 
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Figure 3: Centrality ranking index of Map Kinase cascade protein species. 
The dashed, horizontal lines represent the mean + standard deviation, mean, 
and mean - standard deviation, respectively. 



We further addressed the question of a centrality rank index which rank 
orders the vertices in a hyperdigraph model of chemical reaction networks. 
This centrality rank index highlighted ten protiens of interest in the Map 
Kinase cascade, five of which are highly utilized in the Map Kinase cascade 
while the remaining five are most likely signaling pathway regulators. 

The authors would like to Professor Hong Qian, Dr. Janet B. Jones- 
Oliveira, Professor Stefan E. Schmidt, and Dr. Thorsten Pfeiffer for their 
wonderful insights and discussions. The figures in this paper where gener- 
ated by Graphviz version 1.13(16) a product of AT&T and Pixelglow Soft- 
ware ( [http://www.pixelglow.com/graphviz ). 
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